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Abstract 

Thermal X-ray emission from cooling flows in ellip- 
tical galaxies indicates that ~ f Mq of hot (T ~ 10 7 
K) interstellar gas cools each year, accumulating ~ 
10 10 Mq over a Hubble time. Paradoxically, optical 
and radio frequency emission from the cooled gas is 
lacking, indicating that less than ~ 10 -3 of the cooled 
gas remains. Many have speculated that the cooled 
gas has formed into relatively invisible low mass stars, 
particularly in the context of massive cooling flows in 
galaxy clusters. We focus here on cooling flows in el- 
liptical galaxies like NGC 4472 where the cooled gas is 
made visible in emission lines from HII regions ionized 
and heated (Tjjii ~ 10 4 K) by stellar ultraviolet ra- 
diation. The low filling factor of HII gas requires that 
the hot gas cools at ~ 10 6 cooling sites within several 
kpc of the galactic center. HII mass slowly increases 
at each site at ~ 10 -6 Mq yr _1 until a neutral core 
develops. Neutral cores are heated (Thi ~ 15 K) and 
ionized (x ~ 10 -6 ) by thermal X-rays from the entire 
interstellar cooling flow. We show that the maximum 
mass of spherical HI cores that become gravitation- 
ally unstable is only ~ 2 Mq . No star can exceed this 
mass and fragmentation of collapsing cores produces 
stars of even lower mass. By this means we establish 
with some confidence that the hypothesis of low mass 
star formation is indeed correct - the IMF is bottom 
heavy, but may be optically luminous. Slightly more 
massive stars 5s 4.5 Mq can form near the effective ra- 
dius (r = 8.57 kpc in NGC 4472) if sufficient masses of 
interstellar gas cool there, producing a luminous pop- 
ulation of intermediate mass stars perhaps with radial 
orbits that may contribute to the stellar H/3 index. 
The degree of ionization in gravitationally collapsing 
cores is sufficiently low to allow magnetic fields to dis- 
connect by ambipolar diffusion. Low mass star forma- 



tion is very efficient, involving ~ 10 6 Mq of galactic 
cold gas at any time, in agreement with observed up- 
per limits on cold gas mass. We discuss the cooling 
region surrounding a typical cooling site and show 
that the total X-ray absorption in cold and cooling 
gas is much less that that indicated by recent X-ray 
observations. Using a mass dropout scheme consis- 
tent with X-ray observations and dynamical mass to 
light ratios, we plot the global H/3 surface brightness 
profile in NGC 4472 and compare it with the smaller 
contribution from HII gas recently ejected from red 
giant stars. The lifetime of cooled gas at each cooling 
site, ~ 10 5 yrs, is too short to permit dust formation 
and perhaps also gas phase formation of molecules. 
Subject headings: galaxies: elliptical and lenticular - stars: 
formation - galaxies: cooling flows - galaxies: interstellar 
medium - X-rays: galaxies 

1. INTRODUCTION AND OVERVIEW 

Strong X-ray emission from luminous elliptical 
galaxies is clear evidence that the hot interstellar gas 
they contain is losing energy. Throughout most of the 
galactic volume, this loss of energy does not result in 
lower temperatures since the gas is continuously re- 
heated by compression in the galactic gravitational 
potential as it slowly moves inward. In this sense the 
galactic "cooling flow" is a misnomer. Ultimately, 
however, in the central regions of the flow the gas 
density becomes large enough for radiative losses to 
overwhelm dynamical compression and the gas cools 
catastrophically. For a typical galactic cooling rate, 
~ 1 Mq per year, the total amount of gas that cools 
in a massive elliptical over a Hubble time is large, 
several 10 10 Mq, a few percent of the total stellar 
mass. 

Remarkably, the amount of cold gas observed in el- 
lipticals, cither in atomic or molecular form, is many 
orders of magnitude less than 10 10 Mq (Bregman, 
Hogg & Roberts 1992). The mass of central black 
holes in bright ellipticals is also relatively small, typi- 
cally less than a few 10 9 Mq (Magorrian et al. 1998), 
so the cooled gas cannot be in the holes. Soft X-ray 
absorption has been observed in some galactic cool- 
ing flows, indicating masses of cold gas comparable to 
the predicted value, but the quantitative significance 
or reality of this absorption is unclear at present. 

In addition to cold gas deposited by cooling flows, 
it is possible that additional cold, dusty gas is oc- 
casionally delivered to the centers of ellipticals as a 
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result of merging with gas-rich companion galaxies. 
While this is plausible for some gas-rich ellipticals 
having dusty clouds or lanes, if merging were an im- 
portant source of cold gas for all massive ellipticals, 
the merging rate would need to be carefully regulated 
in order to maintain the small amount of cold gas 
observed in normal ellipticals. 

For many years the standard theoretical explana- 
tion for this shortage of cooled gas is that it has been 
consumed in forming low mass stars (e.g. Fabian, 
Nulsen & Canizares 1982; Thomas 1986; Cowie & 
Binney 1988; Ferland, Fabian, & Johnstone 1994). 
Such young stars must have low masses since neither 
luminous OB stars nor Type II supernovae have been 
observed in normal ellipticals. Explaining the disap- 
pearance of cooled gas by invoking the poorly under- 
stood physics of star formation may seem contrived 
and the low stellar mass hypothesis has led to some 
ridicule. 

The fate of cooled gas in cooling flows associated 
with clusters of galaxies has received most of the ob- 
servational and theoretical attention because of the 
spectacularly large inferred mass deposition rates, 
M > 100 M Q yr- 1 (Fabian 1994). In addition to 
low mass stars, the apparent soft X-ray intrinsic ab- 
sorption of N ~ 10 21 cm -2 indicates that ~ 10 10 Mq 
of cold gas lies within ^100 kpc of the cluster cores. 
Although enormous, this amount of gas would still 
be only a few percent of the total cooled gas based 
on the estimated M, so low mass stars are still the 
preferred endstate for most of the cooled cluster gas 
(Allen & Fabian 1997). However, as with elliptical 
galaxies, this amount of cold gas should in general be 
detectable in HI or CO emission but has not (O'Dea 
et al. 1994; Ferland, Fabian & Johnstone 1994; Vogt 
& Donahue 1995; Puy, Grenacher, & Jetzer 1999), 
amounting to a colossal discrepancy between expec- 
tation and reality. 

In the discussion that follows we revisit the prob- 
lem of cooled gas from the perspective of the galactic 
cooling flow in NGC 4472, a large, well-observed el- 
liptical galaxy. For relatively nearby ellipticals the 
threshold for radio detection is much lower and the 
ratio of observed to predicted cold gas masses is sim- 
ilar to that of more distant cluster cooling flows; e.g. 
H2 and HI are undetected in NGC 4472 with an up- 
per limit 10 7 M Q (Bre gman, Roberts & Giovanclli 
1988; Braine, Henkel & Wiklind 1988), far below the 
~ 10 10 M© expected. Nevertheless, the intrinsic soft 
X-ray absorbing column in NGC 4472, 3 x 10 21 cm~ 2 , 



and its relatively large covering factor indicates cool 
gas masses far in excess of the radio upper limit. The 
cooling flow in NGC 4472 clearly suffers from a minia- 
ture version of the same cooling flow problems of dis- 
tant cluster flows. But because of its proximity, there 
is a large available body of additional observational 
information for NGC 4472 that make it a more appro- 
priate venue to resolve or constrain theoretical possi- 
bilities for the fate of cooled gas. 

But the main advantage afforded by large, rela- 
tively nearby ellipticals emphasized here is that the 
cooled gas is heated, ionized and therefore illuminated 
by ultraviolet radiation from highly evolved galactic 
stars. We argue that the diffuse optical line emission 
from HII gas at T ~ 10 4 K distributed across the cen- 
tral regions of most or all bright ellipticals is a direct 
tracer of cooled gas deposited by the hot interstellar 
gas. A simple analysis of this HII gas leads to the con- 
clusion that hot phase gas is not cooling into a single 
large cloud of neutral gas, but is cooling at a large 
number (~ 10 6 ) of cooling sites located throughout 
the central regions of NGC 4472. The HII gas pro- 
vides direct observational support for a distributed 
mass dropout that has been assumed by many authors 
in the past based on their interpretation of X-ray sur- 
face brightness profiles (e.g. Thomas 1986; Sarazin 
& Ashe 1989). As the mass of HII increases at each 
cooling site, a neutral core of very low temperature 
(T w 10 - 20 K) eventually develops. We show that 
these neutral cores are very weakly ionized and can 
undergo gravitational collapse even in the presence 
of maximum strength magnetic fields. Evidently this 
collapse results in local star formation. 

Another important conclusion from our study of 
the ionized gas in NGC 4472 is that only a very small 
amount, ~ 1 Mq, of neutral (or molecular) gas can 
accumulate at each cooling site before it undergoes 
gravitational collapse. The small mass of collaps- 
ing neutral cores is an essential requirement for low 
mass star formation. Previous studies (e.g. Ferland, 
Fabian & Johnstone 1994) have shown that the gas 
temperature and Jeans mass are small deep within HI 
or H 2 gas irradiated by X-rays in cluster cooling flows, 
but this does not guarantee that massive stars can- 
not form. For example, the Jeans mass is often very 
low in Galactic molecular clouds but these clouds are 
also the birthsites for massive OB stars. The limited 
mass of cold gas at each cooling site in NGC 4472 and 
other similar ellipticals naturally prohibits stars more 
massive than about 1 Mq from forming. 
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The low mass star formation process we propose is 
also very efficient: the total mass of cold gas at all 
cooling sites in NGC 4472 at any time is very small, 
consistent with observed upper limit of cold gas (< 
10 7 Mo). NGC 4472 is an excellent galaxy for study- 
ing this unique star formation process since very lit- 
tle alien gas or stars have been recently accreted into 
NGC 4472 by a merging process. Dusty, and there- 
fore accreted, gas is confined to within r < 0.05r e 
(van Dokkum & Franx 1995; Ferrari et al. 1999) 

A final advantage of studying cooled gas in ellipti- 
cals like NGC 4472 is that the neutral gas formed in 
the cores of HII regions with temperatures T ~ 10 K, 
only lasts for a time, < 10 5 years, that is too short 
for dust (and possibly many molecules) to form. Al- 
though dust and molecules are not required for low 
mass star formation to proceed, these components 
have complicated previous discussions of cooling flows 
in clusters of galaxies where, we assume, the cooling 
process resembles that in NGC 4472. 

One shortcoming of our presentation - as with 
those of previous authors - is that we cannot reconcile 
the observed soft X-ray absorption in NGC 4472 with 
the small amount of cold gas indicated by null obser- 
vations of HI and CO emission. We assume, without 
much justification, that these contradictory observa- 
tions will be resolved in favor of the radio observations 
and that the soft X-ray absorption can be interpreted 
in another way. 

HII gas in elliptical galaxies can also arise from 
stellar mass loss which is ionized by hot central stars 
(planetary nebulae) and galactic UV radiation. We 
begin our discussion below with an argument that 
cooling flow dropout, not stellar mass loss, is likely to 
be the main contributor to internally-produced HII 
gas mass observed in nearby ellipticals. Then we dis- 
cuss an elementary model for the HII gas at a typi- 
cal interstellar cooling site and infer from this a low 
global filling factor for HII gas within the central re- 
gion of NGC 4472. Next we model the cooling of 
hot gas from ~ 10 7 to ~ 10 4 K with a subsonic flow 
in pressure equilibrium - this flow is useful in esti- 
mating the possible contribution of cooling gas to the 
X-ray absorption. This is followed by a discussion 
of the temperature, ionization level and gravitational 
instability of neutral cores at the centers of HII cool- 
ing site clouds. This part of our presentation fol- 
lows rather closely several previous discussions, but 
serves to illustrate that the more spatially concen- 
trated radiative transfer in spherical geometry still 



allows low temperatures and low mass star formation 
within these cores. We then show from observational 
considerations that magnetic fields play little or no 
role in inhibiting the compression of interstellar gas 
as it cools from 10 7 to 10 4 K and from theoretical con- 
siderations that even the strongest observationally al- 
lowed magnetic fields are unlikely to inhibit the final 
collapse of neutral cores to stellar densities. At the 
end of our presentation we discuss the effects of galac- 
tic gravitational forces and stellar collisions on cooling 
site clouds. Finally, to stimulate further observations 
of optical emission lines, we present a surface bright- 
ness map of NGC 4472 showing all the major com- 
ponents: stars, X-ray emitting gas, and HII gas from 
interstellar dropout and stellar ejecta. 

2. NGC 4472: A PROTOTYPICAL ELLIP- 
TICAL GALAXY 

For quantitative estimates in the following discus- 
sion, we use a specific galaxy, NGC 4472, a massive E2 
galaxy associated with the Virgo cluster. NGC 4472 
has been extensively observed at X-ray frequencies 
with Einstein HRI (Trinchicri, Fabbiano, & Canizares 
1986) and with ROSAT HRI and PSPC (Irwin & 
Sarazin 1996). The radial variations of hot gas density 
and temperature based on these X-ray data are illus- 
trated in Brighcnti & Mathews (1997a). Although 
the outer region of the X-ray image of NGC 4472 is 
distorted, possibly by ram pressure interaction with 
ambient Virgo gas, the azimuthally averaged radial 
variation of electron density in NGC 4472 is typical of 
other bright ellipticals (Mathews & Brighcnti 1998). 

The most likely region for low mass star forma- 
tion in NGC 4472 is the volume within ~ 0.1r e where 
r e = 8.57 kpc is the effective or half-light radius at a 
distance of 17 Mpc. The gas that cools in NGC 4472 
cannot all collect at the origin, nor is it likely that 
most of the cooling occurs at very large galactic radii 
where the radiative cooling rate (oc n 2 ) is inefficient. 
Brighcnti & Mathews (1999a) have shown that if all 
the cooled gas accumulates at or near the very cen- 
ter of the galaxy, r < 100 pc, the remaining uncoolcd 
hot gas there is locally compressed and becomes very 
hot, but this is not observed. Alternatively, if most 
of the cooling and low mass star formation occurs in 
0.1 ^ r/r e < 1, then the extremely close agreement 
between the total mass inferred from X-ray data and 
the known stellar mass in this region would be up- 
set (Brighenti & Mathews 1997a). Finally, there is 
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good evidence from observed gas abundance and tem- 
perature gradients that hot interstellar gas is flowing 
inward within ~ 3r e through the optically bright re- 
gions of NGC 4472 (Brighenti & Mathews 1999a), so 
it is unlikely that a significant number of low mass 
stars could form at r > 3r e . It is most interesting 
therefore that HII optical line emission in Ha + [Nil] 
lines is observed just in the region of NGC 4472 where 
low mass star formation is most expected, r < 0.24r e 
(Macchetto et al. 1996). 

3. SEVERAL SOURCES OF HII GAS 

In addition to interstellar gas cooling from the hot 
phase, cold gas is continuously expelled from stars 
throughout the galaxy as a result of normal stel- 
lar evolution. The total rate that mass is supplied 
by a population of old stars in NGC 4472 is M = 
a*(t n )M* t w 1 M© yr -1 where M* t = 7.26 x 10 11 
M Q is the stellar mass in NGC 4472 and a*(t n ) « 
1.7 x 1CT 12 yr _1 is the specific mass loss rate from 
a single burst of star formation after t n — 13 Gyrs 
(Mathews 1989). The supply of gas from stars is com- 
parable to the rate that gas is observed to cool from 
the hot phase: M = (2^m p /5fcT) L Xybo i w 2.5M 
yr _1 , where m p is the proton mass and L x w 7.2 x 10 41 
ergs s^ 1 is the bolomctric X-ray luminosity of NGC 
4472 at a distance of 17 Mpc. 

Several lines of evidence suggest that most of the 
gas lost from stars in ellipticals is dissipatively and 
conductively heated and rapidly merges with the gen- 
eral hot interstellar cooling flow. Gas lost from orbit- 
ing stars inherits stellar velocities which, when dis- 
sipated in shocks or by thermal conductivity, equili- 
brates to the virial temperature of the stellar system, 
T ~ 1 keV. However, the stellar virial temperature is 
about 30 percent lower than that of the more exten- 
sive dark halo. As cooling flow gas slowly flows in- 
ward from the halo into the stellar region, it is cooled 
by ~ 0.3 keV as it mixes with slightly cooler virial- 
ized gas ejected from local stars (Brighenti & Math- 
ews 1999a). This produces the positive temperature 
gradients observed within a few r e . In addition, the 
iron, silicon and other elements supplied by the stars 
increases the metal abundance in the hot interstellar 
gas as it slowly flows toward the galactic center within 
~ r e , producing negative abundance gradients (Mat- 
sushita 1997). These observations indicate that most 
or all of the gas ejected by stars merges with the hot 
gas phase. For simplicity, in the following discussion 



we ignore the HII contribution from stellar mass loss, 
but return to this question in §11. This is contrary 
to the hypothesis of Thomas (1986) that gas ejected 
from stars remains largely neutral and collapses into 
(very) low mass stars without joining the hot phase. 

The assimilation of stellar ejecta into the hot in- 
terstellar gas is greatly accelerated by dynamical and 
thermal processes resulting from the orbital motion of 
mass-losing stars through the cooling flow (Mathews 
1990). Raylcigh- Taylor and other instabilities shred 
the ejected gas into many tiny cloudlets, greatly in- 
creasing the surface area presented to the hot cool- 
ing flow gas and their rapid dissipation by conductive 
heating. In addition, neutral clumps of gas expelled 
from stars always have ionized outer layers which are 
easily ablated and reformed; this results in a rapid 
and complete disruptive heating of the clump. In con- 
trast, cold gas produced as gas cools directly from the 
hot interstellar phase is necessarily formed in the local 
rest frame of the cooling flow gas so the violent dy- 
namical instabilities that accompany stellar mass loss 
are not expected. After ~ 10 6 years, however, the 
denser cooling region may begin to fall in the galac- 
tic gravitational field (see §10), possibly leading to 
some disruption at the cloud boundary (Malagoli et 
al. 1990; Hattori & Habe 1990). Assuming that ra- 
diative cooling from the hot interstellar phase occurs, 
as gas cools through HII temperatures it is thermally 
protected by surrounding gas at intermediate tem- 
peratures (10 4 < T < 10 7 K) where the thermal con- 
ductivity is very low. The global kinematics of the 
two HII gas components of internal origin are quite 
different. HII regions produced by stellar ejecta will 
initially tend to mimic local random and systematic 
stellar motions while HII gas arising from cooling gas 
will initially share the velocity of local hot gas. 

A third source of HII gas in ellipticals are the ion- 
ized parts of gas acquired in recent merging events 
such as the small dusty clouds within ~ 0.05r e of the 
center of NGC 4472 (van Dokkum & Franx 1995). 
This gas is spatially disorganized and is dynamically 
unrclaxcd. Dust is another clue of its external ori- 
gin since gas formed by cooling from the hot phase 
should be nearly dust-free due to sputtering (Draine 
& Salpeter 1979: Tsai & Mathews 1995; 1996) and 
may not have time to grow dust in the gas phase (§7). 

Our interest here is with the HII component pro- 
duced as gas cools from the hot phase and we assume 
that this component dominates the optical line emis- 
sion in NGC 4472. 
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4. THE INVERSE HII REGION 

A small cloud of HII gas that has cooled from the 
hot interstellar medium is photoionized by stellar UV 
radiation arriving at its outer boundary; this is the 
inverse geometry of normal HII regions ionized by a 
central star. We suppose that the HII cloud is spheri- 
cal and that the electron density n e and temperature 
T = 10 4 K are uniform throughout the HII gas. The 
spatial uniformity of the HII density is essentially un- 
affected by small local gravitational fields due to inter- 
nal stars, the neutral core in the cloud if one exists, 
or the HII gas itself. The mass of these HII clouds 
located at the centers of local cooling sites slowly in- 
creases with time. The first step in understanding the 
evolution of HII clouds is to determine the maximum 
size and mass that can be ionized by stellar UV in 
the central regions of NGC 4472. This size depends 
on the HII electron density and the mean intensity of 
galactic UV starlight J uv (r) . 

The intensity of ionizing radiation can be deter- 
mined by an appropriate integral over the galactic 
stellar distribution. For this we assume a de Vau- 
couleurs stellar distribution similar to that in NGC 
4472, with an effective radius of r e = 8.57 kpc and 
an outer maximum stellar radius of 25r e . The stellar 
density and mass are given to a good approximation 

by 

P* = p (6 4 r/r e )-°- 855 exph(6 4 rAe) 1/4 ] 

and 

M(r) =M o7 (8.58,[6 4 r/r e ] 1 / 4 ) 

where b — 7.66925 and 7(0, z) is the incomplete 
gamma function (Mellier & Mathez 1987). 

If the de Vaucoulcurs distribution extends to in- 
finity, the total mass would be M t = M o r(8.58) = 
1.6582 x 10 4 M o where M a = 16tt Po (r e / '6 4 ) 3 . It is nat- 
ural to normalize the density coefficient p Q to fit the 
de Vaucouleurs distribution for NGC 4472 in the re- 
gion 0.1 < r/r e < 1 where the X-ray and stellar mass 
determinations agree, i.e. p Q — 3.80 x 10~ 18 gm/cm 3 . 
When the stellar distribution is truncated at 25r e the 
total mass 6.97 x 10 7 M is only about 1 percent less 
than an infinite stellar distribution having the same 

Po- 

At every radius 

x = r/r e 

in the de Vaucoulcurs distribution the mean stellar 
column density J can be found by integrating over 



solid angle, 




Here p = cos# and £ is the distance back along a 
ray at angle 9 extending through the de Vaucouleurs 
stellar distribution p*(x). The mean column density 
J(x) (in Mq of stellar matter cm -2 ster" 1 ) is calcu- 
lated using the exact solution for the mass intensity 
along each ray and averaging over solid angle with 
192 Gaussian divisions. Then J(x) is converted into 
an approximate mean intensity of ionizing photons 
J(x) by using the conversion factor 7.3 x 10 40 ion- 
izing photons per second per solar mass of old stars 
(Binette et al. 1994; Brighcnti & Mathews 1997b): 

J(x) = 7.3 x 10 40 J(x) cm- 2 ster" 1 s" 1 . 

The local density of ionizing photons is then rn v h {%) = 
4irJ(x)/c. A plot of p*(x) and J(x) in NGC 4472 is 
shown in Figure 1; the variation of J{x) (dashed line) 
and therefore rii P / l (x) is much flatter than the stellar 
density p*(x) (solid line). 

The electron density in HII clouds n can be found 
from 

n(x) = (nT) ism /T HII 

where Thii = 10 4 K and (nT)i sm is the pressure of 
the hot ISM that reproduces the X-ray observations 
of NGC 4472. Plots of n ism (x) and T ism (x) can be 
found in Brighenti & Mathews (1997b). 

The Stromgrcn condition for the maximum radius 
r s of a spherical, fully ionized HII cloud is approxi- 
mately 

C = n 2 a B (4/3)7rr 3 (1) 

where C is the total number of ionizing photons inci- 
dent on (and absorbed by) the cloud per second and 
a B = 2.6 x 1(T 13 (T/10 4 )- - 8 cm 3 s" 1 is the Case B 
recombination coefficient. For an isotropic radiation 
field, a good approximation near the galactic center 
C = AirrlF, where T = 7r J is the flux of ionizing pho- 
tons incident on the surface of the HII cloud. This is 
equivalent to C = {I / A)n ip hcATtr 2 s which is reasonable. 
The Stromgren radius is then 

3ttJ(x) 
r * = — 2 • 

The density of HII gas is p = nMf p where f p = 
5p/(2 + p) = 1.20, assuming p = 0.61 for the molec- 
ular weight. The total mass of the Stromgren sphere 
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-irr 3 s nm p f p = 



36Tr 4 J 3 m p f p 



n- J a 



where m p is the proton mass. Some imprecision is ex- 
pected since we have ignored those ionizing photons 
that pass through the HII cloud unabsorbed. How- 
ever, calculations of the transfer of ionizing radiation 
in the inverse HII region indicate that Equation (1) 
is accurate to <, 5 percent. 

Figure 2 illustrates the radial variation of electron 
density n e = (p/m p )(2 + /x)/5/x in HII gas (solid line) 
with galactic radius in NGC 4472 and the correspond- 
ing local inverse Stromgrcn radius r s (long dashed 
line). Within the radius where Ha is observed in 
NGC 4472, x = r/r e < 0.24 we find r s w 0.3 - 0.8 
pc, n e « 20 — 90 cm~ 3 , and the mass of a typical 
Stromgren cloud is m s w 2 Mq. The radial col- 
umn density in an HII Stromgren cloud is typically 
N s = n e r s w 1.2 X 10 20 cm" 2 . 

Hot gas is assumed to be cooling at numerous sites 
throughout this central region of NGC 4472 and the 
cooling is made visible by optical line emission from 
the HII regions. The mass of any particular HII cloud 
increases slowly with time, supplied by local cool- 
ing from the hot gas phase or by dissipative merging 
of clouds. Presumably, new HII clouds are contin- 
uously forming from the cooling interstellar gas at 
newly-formed cooling sites and old sites and associ- 
ated clouds are disappearing. But we suppose that 
the mean age of cooling sites is long compared to 
the time required for typical HII clouds to reach the 
Stromgren mass; in this case the average cloud can 
be approximated with Stromgren parameters. Notice 
also that r s < r so that even the largest HII clouds 
are very small compared to their distance to the cen- 
ter of the galaxy. 

5. GEOMETRY OF HII AND COLD GAS 

We propose that most of the extended Balmer line 
emission in ellipticals arises from a multitude of HII 
clouds at or near their Stromgrcn radii. If so, the to- 
tal volume within clouds occupies only a tiny fraction 
fp of the galactic volume within the Ha-emitting re- 
gion of NGC 4472, r < 2 kpc. The filling factor f F 
can be estimated by comparing the total volume of 
HII required to produce the observed Balmer line lu- 
minosity to the apparent volume from which optical 
line emission is observed. 

In most optical observations, such as those of Mac- 



chetto et al (1996), Ha (6562 A) is blended with 
two nearby [Nil] lines at 6584 and 6548 A. The 
total flux observed by Macchetto et al. (1996) in 
all three lines is F Unes = 17.30 x 10~ 14 ergs cm~ 2 
s _1 . Observations at higher resolution reveal that the 
F([NII] 6584)/F(Ha) w 1.38 and F(6584)/F(6548) = 
3.0. Combining all these ratios, and adopting Case 
B conditions F Ha /F H p — 2.86, the H/3 flux from 
NGC 4472 is F Hf3 = 2.13 x 10~ 14 ergs cur 2 s" 1 
and the total luminosity from all HII emission is 
L H f3 = <±itD 2 F H p = 7.34 x 10 38 erg s~\ assuming 
a distance of D = 17 Mpc to Virgo. 

How many dust-free Stromgrcn clouds arc required 
to produce this total luminosity? The R(3 luminosity 
of a single Stromgren cloud is 



^ s = n 2 e^(47r/3)r 3 



1.5 x 10 31 n 2 rf pc ergs s" 



where ep = 1.0 x 10~ 25 erg cm 3 s^ 1 is the H/3 emis- 
sivity at T = 10 4 K. For typical values of n e and 
r spc (in parsecs) in the central galaxy x < 0.25, 
£p w 3 - 7.5 x 10 33 ergs s^ 1 . Therefore, about Af cl = 
10 5 — 10 6 Stromgren clouds are required to account 
for the Balmer line luminosity observed. The HII fill- 
ing factor is found by comparing the volume of all HII 
gas V c i = Lhp/ (n e ) 2 eHf3 = 2.5 x 10 60 cm 3 (assuming 
(n e ) = 50 cm~ 3 ) with the total volume of the H/3- 
emitting region, V tot = (4/3)7r(0.24r e ) 3 = 1.1 x 10 66 
cur 3 . The filling factor of HII gas / F = 2x 10~ 6 is 
very small, consistent with our proposition that HII 
emission arises from many small clouds and with ear- 
lier estimates of fp (Baum 1992). If ~ 1 M & of hot 
gas cools each year in NGC 4472, then the mass of 
each cloud will grow quite slowly, ~ 10~ 6 — 10~ 5 M Q 
yr _1 , requiring ( s ~ 2 x 10 5 — 2 x 10 6 years to form 
a typical Stromgren cloud. The total mass of all the 
HII emitting gas in NGC 4472 is M u = {n)Mf F V c i « 
1.2 x 10 5 M Q , similar to values in the literature but 
here evaluated using a consistent physical model. 

A small filling factor also implies that each HII 
cloud is exposed to the unabsorbed stellar UV emis- 
sion from the entire galaxy, provided the cloud sys- 
tem is approximately spherical. The "optical depth" 
for intersecting a Stromgren cloud across the opti- 
cal line-emitting region within r t — 0.24r e is r = 
TTr 2 r t Af c i/V t ot « 0.006 - 0.06. Since t <C 1 the clouds 
do not shadow each other. In reality r could be larger 
(i) if the typical cloud crossection is much less than 
7rr 2 (r cx r^ 1 ) or (ii) if the cloud system were not 
spherical; a disk-like configuration could result from 
galactic rotation. In any case, the assumption that 
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each HII cloud is exposed to the full, unabsorbed stel- 
lar UV emission is likely to be a reasonably good ap- 
proximation. 

Combining previous results, the average column 
depth that HII gas presents to X-radiation through- 
out the galactic core, N <~ N s t ~ 10 18 — 10 19 cm~ 2 , is 
much less than the value N <~ 3 x 10 21 cm -2 that best 
fits the observed X-ray continuum (Buote 1999). The 
size that an HII cloud presents to absorbing X-rays 
is larger than r s since we have ignored the extended 
cooling region around each cloud with temperatures 
between 10 6 and 10 4 K in which X-rays can still be 
absorbed. This assumption will be justified below. 

The total mass of HI or H2 gas observed in the cen- 
tral regions of NGC 4472, M coU < 10 7 M Q , is another 
potential source of X-ray absorption. If this mass of 
cold gas were arranged in a disk of thickness of the X- 
ray absorbing column N = 3 x 10 21 cm~ 2 , located in 
the galactic core and oriented with its symmetry axis 
along the line of sight, it would have a radius < 370 
pc, somewhat larger than the faint patch of dust ob- 
served by van Dokkum & Franx (1996). However a 
cloud of size 370 pc obscures only ~ 0.007 of the total 
X-ray luminosity of NGC 4472 and would therefore 
produce negligible X-ray absorption. The true X-ray 
absorption is very probably much less than 3 x 10 21 
cm~ 2 . The observation of NGC 4472 by Buote us- 
ing the ~ 4' beam of ASCA also included the nearby 
gas-rich dwarf irregular galaxy UGC 7636 (Irwin & 
Sarazin 1996; Irwin, Frayer & Sarazin 1997) which 
may be the source of the X-ray absorbing column at- 
tributed to NGC 4472 if its covering factor is suffi- 
ciently large. 

Although it seems likely that interstellar magnetic 
fields are important in the centers of ellipticals (Math- 
ews & Brighenti 1997; Godon, Soker & White 1998), 
it is remarkable that the observed optical line emis- 
sion does not indicate strong fields in the HII gas. 
Typical HII densities in bright ellipticals determined 
from [SII] 6716/6731 line ratios are - 100-200 cm~ 3 
(Hcckman ct al. 1989; Donahue & Voit 1997), similar 
to (or even a bit larger than) the values found here for 
NGC 4472 (Figure 2). (Unfortunately, we have been 
unable to find a determination of the HII density spe- 
cific to NGC 4472.) This suggests that the HII gas 
density is not being diluted by magnetic support, i.e. 
B 2 /8n < 2nkT or B < 70/iG in the HII gas. HII 
densities of ~ 100 are also supported by comparing 
the ionization parameter U = rn p h/n e for pressure 
equilibrium HII gas in NGC 4472 (short dashed line 



in Figure 2) with values that characterize the entire 
observed line spectrum. Within ~ r e in NGC 4472 
\ogU w —3.3, very similar to values of U required to 
reproduce LINER type spectrum typically observed 
in ellipticals (e.g. Johnstone & Fabian 1988); this 
provides an independent check on our HII gas density 
and J near the center of NGC 4472. 

The apparent absence of magnetic support in the 
HII gas is interesting since the hot phase gas is re- 
quired to have fields of at least several /xG at large 
galactic radii to explain Faraday depolarization of ra- 
dio sources and distant quasars (Garrington & Con- 
way 1991). Interstellar fields > 1/iG can be generated 
in a natural way by stellar seed fields and turbulent 
dynamo action in the hot gas (Mathews & Brighenti 
1997). As the gas density increases by ~ 1000 when 
it cools from the hot phase to HII temperatures, a 
field of 1 \iQ would grow to 100 [iQ> if flux is con- 
served, B cx p 2 / 3 . The initial field in the hot gas 
in r < 0.24r e would need to be surprisingly small, 
< 0.7^iG, to evolve into the rather small fields allowed 
in HII clouds, B < 70^iG, implied by typical electron 
densities. Small HII fields can be understood if local 
cooling sites form in interstellar regions having lower 
than average fields; in pressure balance, lower fields 
require higher hot gas densities which cool preferen- 
tially. Alternatively, it is possible that field reconnec- 
tion has been very efficient during cooling, implying 
a disorganized field and considerable stirring motion 
during the cooling process. 

6. COOLING SITE GAS DYNAMICS 

Cooling sites in the hot interstellar gas are initiated 
in regions of low entropy (i.e. low temperature, high 
density) which cool preferentially by radiative losses. 
Entropy fluctuations can be generated by a variety 
of complex events: stellar mass loss, occasional Type 
la supernovae, sporadic mergers with other nearby 
(dwarf) galaxies, and differential SNII heating events 
and outflows that occurred in pregalactic condensa- 
tions. Due to the complicated nature of these inter- 
actions, it is difficult to predict the amplitudes and 
mass scales of the entropy inhomogencities so the de- 
tailed nature of the initial cooling process remains 
unclear. However, once cooling commences, the gas 
flow toward the cooling site may evolve toward a sim- 
ple profile provided entropy fluctuations in the hot 
gas are not too severe over the flow region. 

We now describe such a model for cooling site 
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flow into an individual HII cloud in which the flow 
is largely subsonic and little influenced by magnetic 
stresses. In the subsonic limit we can assume that the 
flow is isobaric. The thermal energy equation is 

de _ Pdp _ P 2 fA 
^ dt p dt m 2 

where P = (k/ pm p )pT is the constant pressure of 
the ambient hot gas, e = SP/2p is the specific ther- 
mal energy, and A(T) is the Raymond-Cox cooling 
coefficient which must be accompanied by the factor 
/ = (2+/z)(4-3/it)/(25jit 2 ) = 0.58. The corresponding 
temperature evolution is 

dT _ Jl p 2 Pf A 
~dt ~ _ 5 k 2 T' 

If the converging flow at a cooling site is regular over 
times 10 6 — 10 7 years, then steady state is a good 
approximation, m = — pvAirr , and 

dr rhk T 

dt 4:irpm p P r 2 

describes the velocity of a fluid element. The spatial 
variation of gas temperature in the cooling site flow 
is found by dividing the equations above: 

dT 8tt p?P 2 m v f Ar 2 = KK[rs f_ 
dr 5 k 3 m T 2 ~ [ >T 2 ' 

We solve this equation by integrating outward from 
the boundary of a Stromgren cloud of typical radius 
r s = 0.3 pc and temperature T = 10 4 K assum- 
ing a uniform pressure typical of the inner galaxy, 
P = 2 x 10~ 10 dynes. Once T(r) is known, the den- 
sity and flow velocity can be found from the remaining 
equations. Cooling site flows for m = 10~ 6 and 10 -5 
Mq yr" 1 are illustrated in Figure 3. The flow is sub- 
sonic at every radius. When A(T) is slowly varying, 
the temperature varies approximately linearly with 
radius, T oc (ifA) 1//3 r. These simple solutions in- 
dicate that the cooling region in which T < 10 6 K 
is at most about twice the Stromgren radius so the 
extended cooling region can only increase the soft X- 
ray crossection of each cloud by 4, not enough to 
account for the indicated soft X-ray absorption. 

We do not consider conductive heating and ther- 
mal disruption of cooling site flows by ambient hot 
gas since this would be incompatible with our ba- 
sic assumption that the hot interstellar gas is cool- 
ing over a large galactic volume. However, once a 



cooling site is formed, it cannot be easily reversed by 
thermal conductivity. The HII cloud that forms at 
the center of a cooling site is surrounded by a pro- 
tective sheath of gas at intermediate temperatures 
in which the thermal conductivity is extremely low 
(k cx T 5 / 2 ) even in the absence of magnetic fields. 
By comparison, gas expelled from evolving red giant 
stars is directly exposed to interstellar gas at T ~ 10 7 
K, greatly enhancing the likelihood of conductive dis- 
sipation of that gas. Finally, thermal conductivity is 
much less important in galactic cooling flows than in 
hotter cluster cooling flows. 

The soft X-ray emission from the cooling region 
shown in Figure 3 is enhanced by the larger density 
but suppressed by the lower mean temperature in the 
flow. For the flows in Figure 3 the X-ray intensity 
in the ROSAT band (0.2 - 2.5 keV) at r = (ignor- 
ing absorption in the HII cloud) is increased by ~ 10 
compared to the same emission from non-cooling in- 
terstellar gas if it occupied the same volume as the 
flow. However, since the volume of the cooling site 
flow is tiny compared to the entire galaxy, the mean 
X-ray intensity incident on the HI core is not sensibly 
increased by emission from the surrounding cooling 
site flow. Nor is external X-radiation incident on the 
HII cloud appreciably absorbed within the cooling site 
flow. 

7. MASS OF NEUTRAL CORES 

As an HII cloud grows in mass at a cooling site, its 
radius must eventually exceed the Stromgren value, 
causing a neutral or molecular core to develop at its 
center. As the cloud and core mass increase further, 
the core will become unstable to gravitational collapse 
if the gas there is sufficiently cold. We shall determine 
the core mass at which collapse first occurs and thus 
determine an upper limit to the mass of stars that 
can form. For simplicity we assume that the core is 
primarily atomic hydrogen. 

Several authors have discussed the possibility and 
implications of gas-phase dust formation and growth 
in cooled gas in cluster cooling flows (Daines, Fabian 
k Thomas 1994; Fabian et al. 1994; Voit & Donahue 
1995). However, dust may not have time to form or 
grow in NGC 4472 because of the insufficient avail- 
able time. For a lower limit to the dust formation 
time consider assembling graphite grains with carbon 
atoms that stick with every collision and ignoring car- 
bon atom depletion from the gas phase. Grains of 
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radius a g = 10~ 5 cm that can absorb optical light 
have mass m g = 47ra 3 /? c /3 where p c = 2 gm cm -3 is 
the density of graphite. The rate that a grain mass 
grows is m g — n c v c ira^m c where n c = 4x 10~ 4 n#/ 
is the density of C atoms that move at mean veloc- 
ity v c w (8fcT/7rl2m J ,) 1 / 2 . For n H = 10 5 cm" 3 and 
T = 10 K (see below) the grain formation time is 
t g = Ap c a g /12m p n c v c « 2 x 10 5 years. Since this 
minimum value for t g is comparable to the lifetime 
t s ~ 10 5 — 10 6 yrs of neutral or molecular cores (see 
§8 below), it is unlikely that much dust can form. In 
addition, molecules cannot form in gas phase interac- 
tions without grains (Herbst 1987). Formation of H2 
via the H ion is also problematic in galactic cooling 
flows because of the high probability of H~ photode- 
tachment by galactic starlight. We shall not consider 
these details here. 

An essential assumption in estimating the mass of 
neutral cores is that the accumulation of gas in cool- 
ing site clouds is slow compared to the free fall time 
in the neutral core. If neutral gas accumulates faster 
than it can collapse under self-gravity, then more mas- 
sive and luminous stars can form, contrary to obser- 
vation. Since > 10 5 — 10 6 HII clouds are collectively 
processing the entire cooling rate of the galaxy, about 
10~ 6 — 10 -5 Mq flows into each cloud per year. Neu- 
tral cores of density n c are in pressure equilibrium 
with the surrounding HII and ambient hot interstellar 
gas, n c T c « (nT)i sm , ignoring factors of order unity 
involving molecular weights. Near the center of NGC 
4472 the hot gas pressure is (nT) ism w 10 6 K cm" 3 
so the gas density in the neutral cores, 

n c » {nT) ism /T c w 10 5 (T c /10 K)" 1 cm" 3 , 

is large, similar to star-forming cores in Galactic 
molecular clouds. If the cores become gravitationally 
unstable, they collapse in td y ~ (3ir/32n c m p G) 1 ' 2 w 
10 5 years. Since the typical cloud takes a longer time 
~ 3 x 10 5 — 3 x 10 6 to form, the required slowness 
seems to be guaranteed, but not by a wide margin. 

As the masses of molecular cores slowly increase, 
they must eventually become gravitationally unstable 
and form into stars. The condition for this to occur 
can be estimated using the virial theorem with surface 
pressure, 

3 kT c 3GM 2 4 , 

2- -Mhi = t — + 3P-^r 3 2 

2p n m p 5 r c 3 

where P « 2 x 10~ 10 dynes is the ambient interstel- 
lar pressure near the center of NGC 4472. Although 



the temperature and density within the cold core are 
expected to develop radial gradients due to gravity 
forces, in this simple virial model the temperature 
and density within the HI core are uniform, but the 
core pressure is enhanced by self-gravitation. Mag- 
netic forces and rotation are ignored. 

Consider a family of equilibrium models for the HI 
core in which the core mass Mm and external pres- 
sure P are held constant but the core radius r c and 
internal temperature T c are allowed to vary. For very 
large cores, the gravity term in the virial equation is 
small and r c oc T c 1//3 ; for very small cores, the surface 
pressure term is negligible and r c oc T c _1 . Therefore a 
minimum temperature T m exists where dT c /dr c = 0. 
If the true temperature of an HI core is less than T m , 
no equilibrium model is possible and the core must 
collapse. 

It is interesting therefore to examine the locus of 
critically stable models with T c — T m for a wide va- 
riety of core masses, holding the external pressure 
fixed. Differentiating the virial equation and setting 
dT c /dr c = provides an expression for the core radius 
when T c = T m : 



(GMh 
\ 20irP 



1/4 



For each T c > T m there are two radii that satisfy virial 
equilibrium, but only the larger of these roots is stable 

1 /3 

and physically relevant. For this solution, r c oc T c 
holds at large T c , but as T c — > T m the radius decreases 
to r m where gravitational collapse commences. For a 
given core mass and external pressure, the density and 
temperature of these marginally stable configurations 
are 

3M HI ( GMjjj^ ' ' 
\ 20ttP 



Pra = 



4tt 



(3) 



and 



where 



(4) 



K = 



(20tt) 



1/4 



+ 



47T 



3(20tt) 3 / 4 



= 0.751. 



Of most interest is the critically stable core mass for 
any temperature T m : 



M 



HI,m 



(kT m /p n m p ) 2 

K 2pl/2 G 3/2 
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If neutral or molecular cores can cool to < 10 K, col- 
lapse will occur and all stars that form must be less 
massive than ~ 1.56 M© - perhaps very much less 
depending on subsequent mass fragmentation. Since 
the mass of each HII cloud increases slowly compared 
to the free fall time of the neutral core, the core mass 
will not increase further as it collapses. Low mass star 
formation is inevitable. 

The HI number density is related to the core den- 
sity by n = (p/m p )(4-/i n )/(3/i n ) where /i„ = 1.27 is 
the HI molecular weight. For critically stable clouds, 
the neutral column density across the core radius de- 
pends only on the external pressure, 

N =n r 3 4-Mn p0^V /2 

47TTO p 3^L n \ G J 

( P \ V2 

= L1 X ~ — — rrrr- — Cm" 2 (6) 

V2 x io- 10 d y y v ; 

and is independent of the core mass Mhi- Fi- 
nally, the pressure in marginally unstable HI cores 
is P m = AP for all Mhi; this core overpressure rela- 
tive to that in the ambient hot interstellar gas is due 
to the additional gravitational compression unique to 
this region. 

Figure 4 illustrates the mass and radii of the HII 
and HI components of marginally stable clouds as a 
function of core temperature T m with external pres- 
sure P = 2 x 10~ 10 dynes similar to that in the inner 
regions of NGC 4472, r 5S 0.1r e . The minimum sta- 
ble core mass varies inversely with molecular weight, 
Mni.m oc Hn 2 - If hydrogen is molecular, /i#2 = 2.33 
and the leading coefficient in Equation (5) decreases 
from 1.56 to 0.44. Masses for critically stable cloud 
configurations having molecular cores are shown in 
Figure 4 with heavy lines. Also shown in Figure 4 is 
the number density inside critically stable HI cores. 

Equation (4) for the minimum temperature for 
stable core configurations also holds at every radius 
within the core, even if the core density is not con- 
stant with radius, provided we consider a series of 
homologous cores of constant mass in which all radii 
scale linearly and gas density scales as cx r~ 3 . 



8. TEMPERATURE AND STABILITY OF 
NEUTRAL CORES 

The gas temperature within the neutral cores is 
determined by a balance of heating by absorption of 
X-rays or starlight and cooling by far infrared line 
emission. Thermal X-rays are produced by hot cool- 
ing flow gas throughout the galaxy. Infrared line emis- 
sion results from thermal excitation of low-lying fine 
structure levels. If the mean temperature in the core 
(T) based on these physical processes is less than the 
minimum temperature for gravitational equilibrium 
T m , then gravitational collapse and star formation 
will occur. In estimating the temperature in HI cores, 
approximations will be used that tend to maximize 
(T) ; this will lead to an overestimate of the critically 
unstable core mass that collapses. The procedure we 
use is similar to that in other recent papers (O'Dea 
ct al. 1994; Ferland, Fabian & Johnstone 1994; Vogt 
& Donahue 1995), but our treatment of the radiative 
transfer is spherical and we include the possibility of 
UV heating and absorption of X-rays in the HII cloud 
envelope. 

8.1. Photoionization and Heating by Stellar 
UV 

Absorption of stellar UV radiation by those ele- 
ments having ionization potentials less than 1 Ryd- 
berg - C, Fe, Mg, Si, etc. - can be important in 
HI heating and ionization in some circumstances. To 
estimate this possible energy source in HI cores, we 
consider only the most abundant contributor to HI 
heating, carbon ionization. Carbon can be ionized 
by radiation in a narrow radiation window, 11.26 - 
13.6 eV, beyond which radiation is absorbed in the 
surrounding HII region. 

The importance of C — > C + as a source of heating 
and ionization in HI cores can be estimated by balanc- 
ing ionization and recombination rates, just as in an 
inverse HII region. For an HI core of total radius rni 
we expect that carbon will be largely ionized in an 
outer region rc+ < r < rui where rc+ is determined 

by 

C c « n e nc+a c {T)(4/3)Tr(r HI - r% + ). 

Here n e w nc+ ~ nc ~ A c nn is the number density 
of (fully) ionized carbon and A c = 3.62 x 10~ 4 is the 
solar carbon abundance by number relative to hydro- 
gen, The C + + e~ — > C recombination coefficient is 
a c = 1.4 x 10- 10 T-°- 607 cm 3 s" 1 (Hollenbach & Mc- 
Kee 1979; 1989). The equation above can be solved 
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for the radius of the carbon ionization front, 
r C + = r H i(l-j) 1/3 

where 

l{r H i,T) = — — ~~~\2 • 

r HI (A c n H ya c 

If 7 > 1 then rc+ = and carbon is ionized through- 
out the HI core. 

The mean photon intensity of radiation in the car- 
bon ionizing window at any galactic radius J c {r) 
(photons cm~ 2 s _1 ster -1 ) can be estimated from the 
observed flux at a wavelength at the center of the C 
ionization window, A c « 1100 A. For NGC 4472 this 
flux is F c = 2.0 x 10~ 12 ergs cm~ 2 s _1 (Brown, Fergu- 
son, & Davidsen 1995; 1996). At a distance of 17 Mpc 
the luminosity in the carbon window is L c = 6.9 x 10 40 
ergs s _1 or 9.5 x 10 28 ergs s _1 per solar mass in stars. 
This last ratio can be used to convert the mean stellar 
mass surface density J(r) (M cm~ 2 stcr -1 ) to the 
mean intensity of carbon window radiation at any ra- 
dius in NGC 4472, J c (r) = hv c J c (r) = 7.3 x 10 40 J(r) 
= 1.30 x 10 _12 J r (r). Using values appropriate to the 
central regions of NGC 4472, r < 0.1r e , 

/ rr \ 2.607 / N -1 

where P = 2 x 10~ 10 dynes is assumed. We show 
below that the fraction of the HI cores in which carbon 
is ionized is very small, 7 <C 1, so heating by stellar 
UV is unimportant. 

8.2. Photoionization and Heating by Ther- 
mal X-rays 

The bolometric X-ray energy density u x (r) and 
mean intensity J x (r) can be found at every radius in 
the galaxy by averaging the intensity of hot gas ther- 
mal radiation over solid angle. The intensity from 
each direction is found by integrating the equation of 
transfer over the X-ray emissivity e x (n e ,T, z) which 
depends on the local hot gas electron density, temper- 
ature and abundance. The bolometric X-ray emissiv- 
ity for solar abundances can be defined as in Suther- 
land & Dopita (1993), 

(ergs cm~ 3 s _1 ster -1 ) where n p is the proton den- 
sity and fi — 0.61 is the molecular weight of fully 



ionized gas. To evaluate the emissivity e x (r) we 
use the radial variation of hot gas density n e (r) and 
temperature T(r) for NGC 4472 from Brighenti & 
Mathews (1997a), assuming solar abundance. The 
mean intensity J x {r) is found by averaging the X- 
ray intensity over solid angle with 192 Gaussian di- 
visions. The mean bolometric X-ray energy density 
u x (r) = A^J x (r)jc in ergs cm~ 3 is plotted in Figure 
1. In the central regions of the galaxy [r < 0.1r e ) we 
find u x « 10~ 14 ergs cm -3 so that u x ~ Q.lu ip h- 

The specific X-ray photon intensity <fi(E) (photons 
cm~ 2 s _1 kcV -1 ) incident on the cloud is determined 
with the MEKAL option of the XSPEC software code. 
We assume that the X-ray spectrum in NGC 4472 can 
be adequately characterized by an isothermal 1 keV 
plasma with solar abundances. The mean intensity 
at photon energy E incident on the cloud located at 
galactic radius r can be normalized to the bolometric 
energy density u x {r) found earlier: 

j(E ) = CUx ^ E ^ E ^ er § s 
[ ,T) 4tt fE(f>(E)dE cm 2 ssterkcV 

The intensity incident on the outer HII bound- 
ary of a cloud is assumed to be isotropic, 1(E) = 
J(E). The mean intensity at the HII-ISM boundary 
is slightly lower than the ambient J(E) due to absorp- 
tion within the cloud. In computing the X-ray trans- 
fer into the cloud, the incident intensity from 0.01 
to 10 keV is divided into 30 logarithmically-spaced 
energy bands. Only a limited energy range is rele- 
vant. Photons of energy < 0.03 keV cannot penetrate 
the HII gas surrounding HI cores while photons of 
energy > 4 keV pass through critically unstable neu- 
tral cores without absorption. As X-radiation pen- 
etrates the cloud it is absorbed with a crossection 
a(E) modeled after that of Morrison & McCammon 
(1983). The neutral cores are divided into twenty 
equally spaced computational zones and an addi- 
tional zone is assigned to the HII gas. The mean X- 
ray intensity J(E,r) and flux F(E,r) is determined 
at every zone boundary and at the center of the cloud 
for each of the 30 energy bands of mean energy E m . 
The equation of transfer within the cloud is solved to 
find the mean intensity J(E m , r k ) and flux F(E m ,r k ) 
using the same 192-order Gaussian divisions routine 
described earlier, but now including absorption. The 
mean photon intensity J(E m ,rk) is shown in Figure 
5 at three radii - thii = 0.206 pc, rni = 0.041 pc, 
and r = - in an HII-HI cloud having core mass 
Mhi — 3.5 Mq. For this particular neutral core only 
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X-rays with energies 0.1 <j E < 4 keV are absorbed 
in the core. 

The total X-ray heating rate within the radial zone 
between r k and r k +i is 



H k = e(x) 



\F k+1 4irr 2 +l - F k ATvr 2 \ ergs 



(47r/3)(r3 



fe+i 



sec cm° 



(8) 



where F k = T,F(E m ,r k )AE m is the integrated flux at 
r k . For simplicity we assume that the energy of Auger 
electrons is equal to the photoionization threshold for 
each element. Only a fraction e(x) of the total X-ray 
energy absorbed is available to heat the gas. Each 
photoionization by a photon of energy E produces an 
energetic primary electron with energy E — \ where 
X -C E is the ionization potential. This primary elec- 
tron subsequently ionizes many additional H atoms, 
creating on average about E/ Ae secondary electrons 
where Ae = 0.037 keV (Shull & van Steenberg 1985; 
Glassgold, Najita, & Igea 1997). Collisional excita- 
tion of atomic levels of H and He by these secondary 
electrons produces line radiation, most of which es- 
capes and does not heat the gas. Only when the elec- 
tron energy drops below the excitation threshold for 
La (10.2 eV) will all its energy be converted to heat, 
usually by Coulomb exchanges among free electrons. 
The heating efficiency e(x) depends on the ionization 
level x — n(H + )/n(H) which is found by balancing 
rates for recombination and ionization: 

x 2 n(H) 2 a B (T) = (1 - x)n(H)( 

where n(H) is the total number density of both 
neutral and ionized hydrogen and as(T) = 2.91 x 
10- n [T(r)/10K]-°- 602 cm 3 s" 1 is the low tempera- 
ture Case B recombination coefficient (Ferland 1993). 
The ionization rate per neutral hydrogen is 



CM = E 



4nJ(E m ,r)AE n 
1.60 x 10- 9 £ m 



■<r{E m ) 



a7 



where 1.60 x 10 9 converts keV to ergs. The heating 
efficiency in Equation (8) 

e{x) = max{0.10, 0.9971[1 - (1 - ^ ^i.sieS]} 

is an approximate fit to the Monte Carlo results of 
Shull & van Steenberg (1985) but with a constant 
value at very low x (Xu & McCray 1991; Voit & Don- 
ahue 1995). 



8.3. Cooling of HI Cores by Line Emission 

Most of the cooling in the HI gas is due to excita- 
tion of fine structure levels: 370 and 610 /im emission 
from CI and 63 and 145 /xm emission from OI (Hollen- 
bach & McKee 1979, 1989; Ferland, Fabian, & John- 
stone 1994). In estimating the cooling rate of the HI 
gas we determine the exact occupation density of the 
three lowest fine structure levels of CI and OI, taking 
HI collisional rate coefficients and A values from Hol- 
lenbach & McKee (1979; 1989). By considering only 
these dominant cooling levels, we underestimate the 
cooling rate which conservatively overestimates the 
HI temperature and the mass that becomes gravita- 
tionally unstable. 

The principal cooling lines of CI and OI become 
optically thick within the HI core. We allow for this 
by determining the mean escape probability e for each 
line. For each pair of levels (u, I) the line center opti- 
cal depth across the core radius for a Doppler broad- 
ened line is 



1.74 x 10~ 18 NX 3 ul A e iA u , 



9u 
9l 



(T) 



1/2 



where N is the column density in the HI core, \ u i 
is the wavelength of the transition in ^m, A u i is the 
Einstein emission rate, A e i and A e are the cosmic 
abundance and atomic number of C or O, g u and gi 
the degeneracies of the upper and lower levels, and 
(T) is the mass-weighted mean temperature over the 
HI core. Typically r ~ 10 — 100 for strong cooling 
lines. The mean escape probability e for each line is 
found by averaging the mean escape probability from 
a uniform sphere 
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over the Doppler width of the line, 
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Here r D is the optical depth at line center, x = (y — 
v q )/Avd is a dimensionless frequency and Ai>d = 
{v /c)(2k(T)IA e m p ) 1 / 2 is the Doppler width. To an 
excellent approximation in the range 10 _ 1 r < 10 6 
the mean escape probability is 

e(-r ) = -,,!■/ \ where /(r G ) 
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which allows e(r ) to approach unity as r — > 0. The 
mean escape probability for cooling line emission is 
only an approximate substitute for a proper radiative 
transfer calculation including temperature gradients. 
Unfortunately, the escape of line emission can be en- 
hanced by turbulent motion which cannot be easily 
estimated. 

8.4. Results: Gravitationally Unstable Cloud 
Cores 

As the mass of the HI core Mhi is varied, a family 
of critically unstable HII-HI cloud models can be gen- 
erated, all with identical ambient parameters similar 
to those in the optical line emitting region in NGC 
4472: external pressure P = 2 x 10 -10 dynes, bolo- 
metric X-ray energy density u x — 1.7 x 1CP 14 ergs 
s _1 and stellar UV J uv = 1 x 10 8 photons cm~ 2 s _1 
ster -1 . The HII temperature is fixed at 10 4 K and 
the radiative transfer is solved ignoring gravity. The 
temperature and density in each HI core model is ini- 
tially set to T m and p. m throughout. As new tem- 
peratures T(r-fc) are found from the energy balance at 
each radius, the density and radius of each HI zone 
is adjusted to maintain uniform pressure AP through- 
out the core, while preserving each zone mass. The 
X-ray radiative transfer is repeated with the new zone 
spacing, producing new temperatures for the next it- 
eration. The core structure converges after 3 or 4 
iterations. 

Figure 6 illustrates the variation of n(r), T(r), 
M(r), N(r), and x(r) within a neutral core of mass 
3.5 Mq. This configuration is just globally unsta- 
ble to gravitational collapse: the mass- weighted mean 
temperature for this core is (T) = 14.6 K while the 
minimum temperature for stable models having this 
core mass is T m = 15.0 K. Table 1 lists core and HII 
properties of several cloud models having different HI 
core masses. For each cloud the global gravitational 
stability is determined by comparing the mean tem- 
perature (T) with the minimum stable temperature 
T m within radius Thi- It is clear from Equation (7) 
that 7 <C 1 for relevant values of Thi and (T) in Table 
1, i.e. heating and ionization by X-rays dominate UV 
photoelectric processes throughout almost all of the 
HI core volume. 

The two temperatures listed in Table 1, (T) and 
T m , refer to global averages over the entire HI core. 
However, these temperatures can be evaluated at any 
interior radius in the core, as discussed earlier. When 
this is done, we find that (T)/T m has a broad mini- 



mum within the HI structure. Figure 7 illustrates the 
internal variation of (T) /T m within HI core mass for 
cores of five total masses. It is seen that HI cores are 
completely stable as long as the total mass does not 
exceed ~ 2.5 Mq. However, when Mhi ~ 3 Mq the 
entire region between about 1.2 and 2.8 Mq becomes 
unstable. Evidently collapse first begins at an inte- 
rior zone of mass ~ 2 M© when the total core mass is 
2.5 < Mhi ^3.0 Mq. Therefore, as a cloud grows in 
mass at at rate ~ 10 -6 M© yr -1 , stars of mass <, 2 
Mq would collapse at the center of a typical cloud ev- 
ery ~ 2 x 10 6 years. If a single star formed it would be 
of type ^FOV on the main sequence. It is also possi- 
ble that collapsing cores fragment into a distribution 
of stellar masses described by some IMF. For exam- 
ple, if the IMF is Salpeter between 0.03 and 2 Mq, 
then the maximum stellar mass that actually forms 
is about 0.95 Mq. The brighter members of such a 
population of stars would be optically luminous. 

The results in Table 1 confirm that HI core tem- 
peratures are indeed as low as the normalization used 
in Equation (5). Because of the spherical concentra- 
tion of X-ray energy in our calculation, the HI gas 
temperatures in Table 1 are larger than those dis- 
cussed by Ferland, Fabian & Johnstone (1994), but 
their conclusions about low mass star formation still 
hold. The maximum mass of stars that can form at 
cooling sites is ^ 2 Mq. But this instability thresh- 
old is clearly overestimated. The unstable core mass 
would be significantly smaller, for example, if some of 
the core gas were molecular, if more radiative cool- 
ing transitions were considered, or if a small amount 
of internal turbulence lowered the optical depths of 
cooling line radiation. 

Is massive star formation also excluded at larger 
galactic radii? To explore this wc performed a se- 
ries of calculations like those in Table 1 but using 
external parameters expected at r = r e = 8.57 kpc 
in NGC 4472: pressure P = 10~ n dynes, X-ray en- 
ergy density u x = 6.1 x 10~ 16 ergs cm -3 and UV 
intensity J uv = 4 x 10 6 photons cm~ 2 s _1 ster -1 . 
With these parameters we find that the lowest grav- 
itationally unstable HI core mass is 4.5 Mq where 
(T) « T m ks 9 K. The lower pressure at r = r e gener- 
ates somewhat larger critically unstable core masses, 
even though the mean HI temperature is lower (Equa- 
tion 5). Fragmentation into a Salpeter IMF (from 
0.03 to 4.5 Mq) during collapse would limit the typ- 
ical maximum stellar mass to ~ 2 Mq (i.e. A5V to 
F0V). Such stars, having highly radial orbits, could 
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contribute to the H/3 absorption line indices observed 
in elliptical galaxies (Gonzalez 1993; Worthey 1994; 
Faber et al. 1995; Worthey, Trager, & Faber 1995) 
which have been interpreted as evidence for a popula- 
tion of young stars of intermediate mass in ellipticals. 
The observed H/3 absorption index usually strength- 
ens toward the centers of ellipticals (in r < 0.5r e ), 
indicating lower mean stellar ages; this is consistent 
with cooling flow dropout which should also be cen- 
trally enhanced. The more prominent H/3 absorption 
seen in ellipticals of lower luminosity may reflect the 
influence of lower interstellar pressure in Equation (5) 
and correspondingly larger unstable core masses. 

9. COMPLICATIONS: MAGNETIC FIELDS 
AND ROTATION 

For a definitive understanding of star formation 
in elliptical galaxies or elsewhere, it is necessary to 
consider the influence of magnetic fields and rotation, 
both of which are intensified during the two major 
compressions that lead to gravitational instability and 
star formation in elliptical galaxies. 

During the initial cooling condensation from T = 
10 7 to 10 4 K, neither rotation nor magnetic energy 
density appear to have strongly restricted the col- 
lapse. Any fragmentation due to rotation during this 
initial cooling phase probably serves only to increase 
the multitude of local cooling sites, each containing 
a small HII region. Remarkably, the frozen-in ampli- 
fication of magnetic energy during this initial com- 
pression has not occurred since optical observations 
indicate that the HII gas pressure equals the pressure 
of the local hot interstellar gas and is therefore undi- 
luted by additional magnetic support. The second 
cooling compression by another factor of <~ 10 3 oc- 
curs as the core gas recombines and the temperature 
falls from 10 4 K to ~ 10 K. During this compression, 
magnetic fields will intensify if they remain attached 
to the weakly ionized plasma and rotation should also 
increase. Rotation during this final stage may lead to 
subsequent fragmentation into even smaller masses, 
although the entire cloud is optically thick in the main 
cooling lines. A helpful interplay between magnetic 
breaking and angular momentum redistribution can 
be expected as stellar densities are approached. In 
this section we discuss a limited aspect of this general 
problem, the loss of magnetic stresses by reconnection 
and ambipolar diffusion as neutral gas recombines and 
compresses to typical HI core densities, ~ 10 5 cm~ 3 . 



Although difficult to prove rigorously, reconnec- 
tion losses could in principle greatly reduce magnetic 
stresses and their inhibiting influence on compression. 
We have already noted that the observed HII gas pres- 
sure in bright ellipticals restricts fields in the HII gas 
to B Sa 70 fiG, suggesting that some reconnection 
losses occurred during compression from local inter- 
stellar to HII densities. Reconnection could be stimu- 
lated by turbulence and disordered fields expected in 
the hot interstellar gas (Mathews & Brighenti 1997) 
which, when compacted to HII densities, result in lo- 
cally anti-directed fields. In the second compression 
phase, as HII gas condenses into HI cores, the density 
increases by another factor of 1000 and the final flux- 
frozen fields would then be < 7000 fiG. Fields this 
large in the HI core would strongly inhibit gravita- 
tional instability and star formation and increase the 
minimum core mass necessary for gravitational col- 
lapse. However, it is possible that reconnection losses 
also dominate this second cooling phase if the field 
is sufficiently disordered (Friaca & Jafelice 1998). If 
reconnection is wonderfully efficient, the field- free es- 
timates we have made here for the maximum unstable 
stellar mass do not require further modification. 

Fortunately, physical conditions in HI cores allow 
for efficient slippage of the magnetic field relative to 
the neutral gas. For densities and field strengths of in- 
terest here, the charged-neutral collision rate is much 
less than the Larmor frequency for both ions and 
electrons, so ambipolar diffusion should prevail over 
Ohmic dissipation or Hall effect (Wardle & Ng 1999). 
As the field moves relative to the (essentially station- 
ary) neutral component in HI cores, ion-neutral drag 
must balance magnetic tension, (av}rnpvd ~ B 2 /Aitr. 
Here a is the momentum transfer cross-section. When 
Vd is less than some transition velocity v tr , charge- 
neutral interactions are dominated by the polarizibil- 
ity a of atomic hydrogen, a ss 2.47r(e 2 a/^ m w) 1/ ' 2 and 
the rate coefficient (av) w 2 x 10~ 9 cm 3 s _1 is in- 
dependent of temperature (Osterbrock 1961; Nakano 
& Umebayashi 1986). When Vd Si v tr the crossection 
is limited by the geometrical size of the particles in- 
volved, a w ita 2 Bl where the Bohr radius as is appro- 
priate for H + or e~ collisions with H (Mouschovias 
& Paleologou 1981). When applied to the equa- 
tion above, the drift velocity Vd is either linear or 
quadratic. For collisions in weakly ionized hydrogen, 
the transition velocity v tr ~ 230 km s _1 is found by 
equating the two crossections. 

The time required for a field to leak through a 
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weakly ionized neutral core t a d ss r/vd in the linear 
regime is 



t 
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where x = n e /n is the degree of ionization. 

To illustrate the efficiency of ambipolar diffusion 
in collapsing HI cores we can exploit the result that 
larger fields diffuse out faster, i.e. t a d oc B~ 2 . We 
begin by considering two simple limits: (1) initially 
uncompressed neutral gas that is supported by mag- 
netic stresses and (2) dense neutral gas after it has 
collapsed to core densities approaching gravitational 
instability. For the initial state, consider a gas which 
has become largely neutral, but is supported by a field 
of maximum strength, B ~ 70/xG, so that its density 
n w 100 cm~ 3 is still that of the HII gas. As long as 
= (B 2 /8ir)/P < 0.5 the column density in such 
a core of mass ~ 1 M Q still resembles that of typi- 
cal HI cores in Figure 6, so x ~ 10~ 6 . Under these 
circumstances the ambipolar diffusion time is 
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This time is much less than the free fall time for field- 
free cores, tff ss 5 x 10 6 yrs, so the strongest possible 
initial fields cannot arrest the initial compression to- 
ward typical HI core densities. 

After the HI gas compresses to the final core den- 
sity but just before gravitational collapse occurs, then 
n w 10 5 cm -3 and a magnetic field ~ lOO^G is com- 
parable with the gas pressure. In this final state be- 
fore core collapse, the ambipolar diffusion time 

(2) 6 ( ^* \ { ^ \ 



yrs 



,3 x io- 7 y Vioo/xGy 

is now longer than the core freefall time, ~ 10 5 yrs or 
the dynamical time ~ 10 6 yrs desired to process cold 
gas rapidly into stars, involving observationally small 
HI and CO masses. Some additional field losses by 
reconnection may be required. 

However, this estimate of the field diffusion time in 
the final core fails to acknowledge that field diffusion 
occurs steadily during the collapse to core densities. 



Suppose for example that the field varies as if it were 
frozen-in, B oc n 2 / 3 , then the ratio of field diffusion to 
free fall time at each density during the collapse varies 
as tad/ tff oc v}l 2 x which increases only modestly as 
n increases by ~ 1000, holding x approximately con- 
stant. (We use the linear drift regime, which generally 
applies.) If this time scale ratio is normalized to the 
initial time ratio above, t a d ^ /tff = 0.0013, then 
t a d ^ /tff = 0.04 when core densities are reached, 
which remains much less than unity. In the absence 
of a detailed calculation, this crude estimate sug- 
gests that magnetic fields can slip away from HI cores 
rapidly compared to the dynamical times involved in 
the gravitational evolution of HI cores. 

When field diffusion is important, however, the 
field is spatially redistributed, not destroyed. In a 
rigorous calculation it would be necessary to follow 
the subsequent evolution of the field which may accu- 
mulate in the outer regions of the core. If left behind 
after star formation, such magnetically enhanced re- 
gions may ultimately become buoyant in the galactic 
field and merge again with the HII gas. In addition, 
the thermal balance in the neutral gas will be altered 
by heating due to field diffusion and reconnection, but 
this heat may be carried away by efficient radiative 
losses. In summary, it appears that magnetic stresses 
may not be an overwhelming problem in restricting 
neutral cores from gravitational collapse. How mag- 
netic and rotation effects influence the subsequent his- 
tory of the star after collapse ensues will not be en- 
tered into here. 



10. DYNAMICAL INFLUENCES ON COOL- 
ING SITE CLOUDS 

Because of their high density relative to the hot 
interstellar gas, HI-HII clouds will tend to fall in the 
gravitational field of the galaxy past the hot gas. As 
an HI-HII cloud falls toward the center of the galaxy, 
the HII envelope should experience a drag force, so 
some motion of the HI core away from the center of 
the HII envelope could be expected. For a conser- 
vative estimate of this, suppose that the HI core is 
accelerated toward the center of the galaxy by the 
galactic gravitational force but the galactic position of 
the HII envelope is held fixed. At a galactic radius of 
r = 0.03r e , the gravitational acceleration toward the 
galactic center is g « GM^(r)/r 2 w 2 x 10~ 6 cm/s 2 
where M t w 7 x 10 9 M at r = 0.03r e . (If the cooling 
flow is rotating, this acceleration will be an overesti- 
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mate.) The time required for the HI core to move a 
distance equal to its own diameter (2r c w 0.08 pc), 
td = {Arc/g) 1 / 2 w 1.5 x 10 4 yrs, exceeds the recombi- 
nation time in the HII gas, t rec — l/(nas) « 1.4 xlO 3 
years. Even in this extreme example which maximizes 
the differential motion of HII and HI, gas can ex- 
change between the HI and HII phases rapidly enough 
to maintain the neutral core near the center of the 
cloud. 

How many galactic stars are momentarily passing 
through HII clouds or their cold cores at any time? 
From Figure 1 the stellar density at r = 0.1r e is 
p* w 2 x 10~ 22 gm cm~ 3 or 3 M Q pc~ 3 . If the average 
stellar mass is 0.3 M , the stellar density is ~ 10 stars 
pc~ 3 . An HI-HII cloud with radius of r HII w 0.21 pc 
(Table 1) has volume ~ 1 pc~ 3 and should contain 
about ten small stars. The time-averaged density of 
stars within the HII region is comparable to that of 
the HII gas itself (hh ~ 2 Mq pc -3 ). Moving at a 
typical orbital velocity, ~ 300 km s _1 , a star should 
pass through the tiny HI core of radius mi = 0.04 
pc only very infrequently, every ~ 10 5 years, caus- 
ing little or no disruption. Finally, we note that the 
current Type la supernova rate required to maintain 
the observed interstellar metallicity, ~ 0.03 SNe per 
10 10 Ls every 100 years, is so low that the galactic 
volume occupied by supernova remnants is negligible. 
Supernovae are not expected to greatly disturb cool- 
ing site clouds. 

11. MASS DROPOUT AND EMISSION LINE 
SURFACE BRIGHTNESS PROFILES 

Optical line emission in ellipticals identifies where 
hot interstellar gas is cooling and forming into low 
mass stars. In this section we estimate the H/3 surface 
brightness in NGC 4472 due to emission from cooled 
hot interstellar gas. We use a dropout procedure that 
is self-consistent, but not necessarily unique. 

One of the least understood but essential compo- 
nents in evolutionary models for the hot interstellar 
gas is the cooling dropout mass profile. The notion of 
distributed cooling mass dropout was first invoked in 
steady state cooling flow calculations to reduce the 
central X-ray surface brightness peak to better fit 
the observations (Thomas 1986; Thomas et al. 1986; 
Sarazin & Ashe 1989; Bertin & Toniazzo 1996). In our 
recent galactic cooling flow models we compare the ef- 
fect of distributed mass dropout on both the apparent 
(observed) temperature and electron density distribu- 



tions. Dropout reduces the apparent temperature and 
increases the apparent density. However, we find that 
it is the dropout contribution to the mass to light ratio 
(dynamically determined from stellar velocities) that 
most strongly constrains cooling dropout mass pro- 
files (Brighenti & Mathews 1999b). One of the sim- 
plest dropout assumptions that adequately satisfies 
these constraints is the recipe employed by Sarazin & 
Ashe (1989); we use it here to illustrate the relation- 
ship between optical line emission and mass dropout. 

Following Sarazin & Ashe (1989), we suppose that 
the hot gas cools at a rate p = qp/tdo that depends 
on the local hot gas density and the cooling time at 
constant pressure, tdo = 5m p kT/2pAp. Here T is the 
hot phase temperature and q is a constant factor that 
characterizes the cooling rate. Assuming that the 
M/Lb of low mass stars is very large, we find that 
q = 1 satisfies reasonably well most of the observa- 
tional constraints for NGC 4472 (Brighenti & Math- 
ews 1999b). We suppose that most of the cooling 
gas is processed through mature cooling sites where 
the typical HII cloud mass is comparable to the lo- 
cal Stromgren mass m s . At each cooling site mass 
is processed into low mass stars at a rate ~ m s /t s 
where t s ~ 10 5 yrs, based on the filling factor ar- 
gument. By continuity, the total rate that cooling 
mass passes through an element of galactic volume 
dV is pdV w dJV(m s /t s ) where dM is the num- 
ber of ^Stromgren clouds in dV. Each Stromgren 
cloud has an H/3 luminosity £p t3 w nepm s / pm p where 
n = (nT)i sm /THii is the electron density of gas at 
Thii = 10 4 K and ep is the Case B H/3 emissivity. 
The mean global H/3 emissivity at radius r due to 
dropout cooling is therefore 

dLp _ dM _ . t s nepm s 
6/3 A-Kr 2 dr dV /3 ' s ^ m s pm p 

= q t s ^Lnl m A(T lsm ) J** (9) 
okl hii cm a s 

Notice that ep is independent of the mean cloud mass 
m s but depends on the time t s for low mass stars to 
form at each cooling site. The H/3 surface brightness 
T,p(R) for this dropout model is found by integrat- 
ing ep along the line of sight through the galaxy at 
projected radius R. 

For comparison, we also estimate the H/3 surface 
brightness distribution from gas recently ejected from 
evolving stars and ionized by galactic ultraviolet ra- 
diation. Ejected envelopes from red giant stars are 
not sufficiently massive to produce clouds of Strom- 
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gren mass. The ejection process is complex and slow, 
lasting ~ 10 5 years. In addition, mass-losing stars 
orbit rapidly through the hot gas. This results in vio- 
lent Rayleigh- Taylor disruption in ~ 10 4 years as the 
ejected gas is decelerated and shocked by the hot in- 
terstellar gas. The huge surface area that the highly 
fragmented stellar ejecta presents to the hot gas pro- 
motes the conductive assimilation of this HII gas into 
the local hot gas on timescale t* con d ~ 10 4 yrs (Math- 
ews 1990). The assimilation time t* con d depends on 
the density of the local interstellar gas, the mean ve- 
locity of mass-losing stars, and other parameters. In 
estimating the H/3 emissivity from stellar ejecta we as- 
sume that these complications can be subsumed into 
the single parameter t* con d which we regard as spa- 
tially uniform. The rate of supply of mass from evolv- 
ing stars is proportional to the local stellar density, 
p* = where a* = 5.4 x 10~ 20 sec -1 is appro- 

priate for a Salpeter type IMF. The local mean H/3 
emissivity due to stellar ejecta is therefore 

dL p ^ _ , ntp erg 

inr-'ar fim p cm 3 s 

where n = {nT)i sm /THii since pressure equilibrium 
is reached on very short timescales. 

In Figure 8 we illustrate the X-ray surface bright- 
ness Y>rosat of NGC 4472 in the ROSAT energy 
band (0.2 - 2.0 keV) computed from the combined 
density and temperature distributions n(r) and T{r) 
determined from Einstein, ROSAT, and ASCA data 
(see Brighenti & Mathews 1997a for details). Using 
these same n(r) and T(r) profiles, we have deter- 
mined the H/3 emissivity from Equation (9) and plot 
the corresponding surface brightness £#/3 (dropout) 
in Figure 8, assuming a Sarazin-Ashe dropout coef- 
ficient q = 1. Y,hp (dropout) is normalized so that 
the total H/3 luminosity within R w 0.25r e matches 
the observed value Lnp = 7.34 x 10 38 erg s _1 from 
this region; with this normalization the total H/3 lu- 
minosity from the entire galaxy is somewhat larger: 
Lh/3 = 11-8 x 10 38 erg s _1 . To achieve this normal- 
ization it is gratifying that the characteristic star for- 
mation time in Equation (9) must be t s = 2.2 x 10 5 /q 
years, very close to the value discussed earlier when 
q = 1. This value of t s normalized to Lu(i{r < 0.25r e ) 
is a consequence of our assumption that most H/3 line 
emission comes from dropout, any additional HII con- 
tribution from stellar mass loss would tend to increase 
t„. 

We also plot in Figure 8 the stellar surface bright- 



ness £* for NGC 4472 using the standard de Vau- 
couleurs relation 

Z*(R/r e ) - -^cxp{-7.67[(i?/r e ) 1 /4 _ 1]} 

which agrees very well with the CCD photometry of 
NGC 4472 by Peletier et al. (1990). The H/3 sur- 
face brightness profile from gas recently ejected 
from stars using Equation (10) is also shown in Fig- 
ure 8 with two assimilation times, i*, ccm d = 10 4 (pre- 
ferred) and 10 5 yrs. These estimates for Y,h/3,*, al- 
though very uncertain, are less than (dropout), 
indicating that most of the H/3 observed in NGC 4472 
is related to mass dropout not stellar mass loss. This 
conclusion is not inconsistent with the observation of 
a few bright, young planetary nebulae in ellipticals 
ionized by central stars (e.g. Arnaboldi, Freeman & 
Gerhard 1998). A considerable amount of observa- 
tional data is available on the radial distribution and 
kinematics of optical line emission in normal ellipti- 
cals (Trinchieri & di Serego Alighieri 1991; Macchetto 
et al. 1996; Zeilinger et al. 1996; Pizzella et al. 1997; 
Trinchieri, Noris & di Serego Alighieri 1997; Goud- 
frooij & Trinchieri 1998; Corsini et al. 1999), but ad- 
ditional observations of optical emission line spectra 
and surface brightness distributions would be useful 
in understanding the physics of mass dropout from 
the hot interstellar gas. 

12. FINAL REMARKS AND CONCLUSIONS 

Models for the evolution of hot interstellar gas in 
ellipticals that are successful in most other respects - 
such as explaining the observed radial distribution of 
hot gas density, temperature and abundance - cannot 
be fully accepted unless they can also account for the 
final disposition of cooled gas. Recent evolutionary 
models for the hot gas in NGC 4472 require that a 
mass M coo i ~4x 10 10 M Q must cool in some fashion 
over the Hubble time (Brighenti & Mathews 1999b). 
There are two serious problems in discharging this 
cold gas from the cooling flow: (1) Although M coo i 
is much less than the total stellar mass M* = 7.26 x 
10 11 Mq, the mass of cooled gas must be deposited 
within the (inner) galaxy in a manner consistent with 
the known stellar mass to light ratio, and (2) The 
cooling process and the conversion of cooled gas into 
low mass stars or some optically dark equivalent must 
be supported with a physically plausible model - this 
is the problem we address here. 
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Our objective in this paper is to further our un- 
derstanding of mass dropout by exploiting features 
unique to elliptical galaxy cooling flows that are gen- 
erally unavailable or less obvious in cluster cooling 
flows. Foremost among the advantages in galactic 
flows is that the the dropout is illuminated by stel- 
lar UV ionizing radiation, appearing as HII emission 
visible throughout the inner galaxy. By this means 
we assert that the cooling must be widespread, dis- 
tributed among many cooling sites. Moreover, the 
cooling rate at each site is sufficiently slow that HII 
clouds there can develop neutral cores that undergo 
gravitational collapse after only ~ 2 M Q of neutral 
gas has accumulated. Only low mass stars can form. 
The conversion of hot gas into low mass stars in galac- 
tic cooling flows is very efficient, with star-formation 
timescales so short, ~ 10 5 years, that dust and prob- 
ably molecules cannot form. Globally, this cooling 
dropout model involves very little mass of neutral gas 
at any time, Af c iouds m s ~ 10 6 M Q and satisfies the 
low upper limits on observed HI + H 2 gas in NGC 
4472 (< 10 7 M ). 

In a similar investigation several years ago, Fcr- 
land, Fabian & Johnstone (1994) studied the thermal 
structure of cold clouds directly illuminated by X- 
rays from cluster cooling flow gas. Specifically, they 
solved the plane parallel transfer of ~ 5 keV thermal 
radiation into a cold, non-magnetic, dust-free slab, 
including a large number of coolants and chemical re- 
actions. They showed that deep within the irradiated 
slab, TV > 10 20 cm -2 , the gas becomes molecular and 
the temperature falls below 10 K where the Jeans 
mass is only a few tenths of a solar mass, implying 
that low mass stars are likely to form. 

Our study here is a complementary extension of 
this idea, including several additional considerations. 
We consider fewer atomic and chemical processes than 
the Ferland, Fabian & Johnston (1994) calculation, 
but many complicated processes involving dust and 
molecules may be irrelevant in galactic cooling flows 
due to the short lifetime of the cooled gas before col- 
lapsing to stellar densities. We have shown that the 
more concentrated spherical radiative transfer of X- 
rays into cooling sites with star-forming neutral cores 
does not increase the core gas temperature appre- 
ciably. Unlike the plane parallel calculation, where 
Jeans unstable regions could in principle merge to 
form more massive stars, in our spherical model the 
total mass of gas in the cold cores prior to gravita- 
tional collapse - < 2 M - absolutely precludes the 



formation of massive stars which are never observed 
in normal ellipticals. 

The simple model proposed here is unlikely to ap- 
ply to every elliptical galaxy. Unlike massive ellipti- 
cals like NGC 4472, the stellar systems in low lumi- 
nosity ellipticals (My < —21) differ in many impor- 
tant details: they are rotationally supported, have 
more nearly isotropic stellar velocities, have dense 
power-law stellar cusps, are disky, etc. (Faber et 
al. 1997). Brighcnti & Mathews (1997b) have shown 
that large, rotationally supported cold gaseous disks 
can form in these low luminosity ellipticals and that 
normal, luminous star formation in the disk may ac- 
count for the systematically higher stellar H/3 indices 
observed in these galaxies (de Jong & Davies 1997, al- 
though their correlation with H/3 has been challenged 
by Rampazzo et al. 1999). Recent HI observations of 
low luminosity ellipticals have detected smooth, dif- 
ferentially rotating disks of HI gas in these galaxies, 
sometimes extending right to the center (Oosterloo, 
Morganti, & Sadler 1998) similar to the cold disks 
predicted by our models. Evidently, star formation in 
large, coherent HI disks in low mass ellipticals differs 
from the model described here in which the HI gas is 
distributed among many isolated, physically distinct 
clouds. 

Small irregular dusty clouds are observed in the 
cores of most (~ 80 percent) bright ellipticals (van 
Dokkum & Franx 1996). Centrally orbiting HII gas 
is also occasionally observed, sometimes rotating in 
a sense opposite to that of the galaxy (e.g. Bertola, 
Buson & Zollinger 1988). In these cases gas has pre- 
sumably been acquired by a merging process unre- 
lated to cooling flow evolution. Since the dusty clouds 
typically have masses in excess of the critical masses 
for gravitational collapse ~ 1 M©, luminous massive 
star formation with a normal IMF can be expected 
as well as Type II supernovae, but at very low rates. 
Normalizing to our Galaxy where there are ~ 10~ 12 
SNII yr _1 per solar mass of molecular interstellar gas, 
the < 10 7 M Q of cold gas in NGC 4472 should gen- 
erate < 10~ 5 SNII yr _1 . Because this SNII rate is 
so small (and their light is masked by starlight in 
the galactic core), it is not surprising that SNII have 
escaped detection in ellipticals. Optical luminosity 
and H/3 emission from normal star formation in dusty 
cores is also deeply buried in starlight. According 
to the data of Tacconi & Young (1986), the disk of 
the Scd galaxy NGC 6946 produces ~ 3.5 x 10 33 ergs 
s _1 of B-band light per solar mass of H 2 gas so the 
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M < 10 7 Mq of cold gas in NGC 4472 should gen- 
erate Lb a < 3.5 x 10 40 erg s _1 . The dusty gas 
cloud in NGC 4472 has a maximum extent of only 
rd ~ 0.91 arcseconds (van Dokkum & Franx 1995) 
which is less than the "break" radius r\, = 2.41" 
in the stellar light. Assuming the B-band surface 
brightness in NGC 4472 is uniform within r^, the to- 
tal projected luminosity produced by galactic stars, 
L B ,proj{rd) ~= Trr^EsiVd) w 4.5 x 10 40 erg s" 1 , ex- 
ceeds LB,d- We conclude that normal stellar forma- 
tion can occur in accreted dusty clouds at the centers 
of bright ellipticals like NGC 4472 without easy de- 
tection. 

The interesting possibility of observing cold gas in 
elliptical galaxies with the 6.4 kcV fluorescent Fe line 
has been discussed by Churazov et al. (1998). How- 
ever, the small star-forming regions discussed here are 
completely transparent to K shell ionizing radiation 
at 7 keV (see Figure 4) so very little 6.4 keV radiation 
would be expected. In any case, galactic cooling flows 
of temperature ~ 1 keV have little emission above 7 
keV. Perhaps the best confirmation of our star for- 
mation model would be to directly observe the young, 
low mass stars and a variety of infrared techniques are 
available for this purpose (Kroupa & Gilmore 1994; 
Joy ct al. 1995; Prestwich et al. 1997). 

Our conclusions can be briefly summarized as fol- 
lows. The final stages of catastrophic cooling in galac- 
tic cooling flows is distributed among a large num- 
ber of cooling sites, each occupying a tiny fraction 
of the galactic volume. The mass of cold gas within 
HII clouds at each site accumulates slowly, allowing 
episodic gravitational collapse of dense neutral cores. 
In the typical interstellar environment of giant ellip- 
tical galaxies, the maximum mass of gravitationally 
unstable regions is only <> 2 Mq (i.e. later than A5V), 
setting an absolute upper limit on the mass of stars 
that can form. Some optical luminosity would be pro- 
duced by these stars even if the collapsing cores frag- 
mented with a normal IMF. Magnetic stresses and 
rotation could complicate this conclusion, but simple 
estimates indicate that magnetic fields diffuse rapidly 
out of collapsing, star-forming cores. This method 
of forming low mass stars involves very little neutral 
or molecular gas and is completely consistent with 
the apparent absence of cold gas in normal ellipticals 
- assuming that the mass of cold gas is more accu- 
rately determined from radio than from X-ray obser- 
vations. If cooling from hot interstellar gas occurs 
at large galactic radii, near r e , then stars of larger 



mass < 2 — 3 Mq can form. Such stars may help pro- 
duce the H/3 absorption features observed in elliptical 
galaxies which have been interpreted as evidence for 
a population of youthful stars of intermediate mass. 
Finally, in elliptical galaxies the star-forming sites are 
illuminated by optical HII line emission so the surface 
brightness of this radiation must be in accord with 
mass dropout models as well as with observations of 
the hot X-ray emitting interstellar gas. 
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Table 1 

GRAVITATIONAL STABILITY OF HI CORES a 
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models. 



h i"Hll is the outer radius of the entire cloud, most of which is ionized. 

c Mass- weighted average temperature in HI core. 

d Minimum temperature for pressure-supported HI cloud cores. 
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Fig. 1. — Radial variation of stellar density 
P*(r/r e ) in NGC 4472 (solid line) and correspond- 
ing mean photon intensity in the H-ionizing contin- 
uum, J(r/r e ) photons cm~ 2 s _1 stcr -1 (dashed line). 
The bolomctric energy density of thermal X-radiation 
from the hot interstellar g clS Ux 

(r/r e ) is shown for 
comparison (long-short dashed line). The radius is 
normalized to the effective radius r e = 1.73 arc min- 
utes or 8.57 kpc at D = 17 Mpc. 

Fig. 2. — Radial variation of HII gas density n e (solid 
line), HII Stromgren radius r s (long dashed line), and 
ionization parameter U = n ip h/n e (short dashed line) 
as functions of galactic radius in NGC 4472. The 
radius is normalized to the effective radius r e = 1.73 
arc minutes or 8.57 kpc at d = 17 Mpc. 

Fig. 3. — Steady state spherical gas flow surrounding 
HII gas of radius 0.3 pc at a typical galactic cooling 
site. Gas temperature, velocity and density are shown 
with solid, long-dashed and short-dashed lines respec- 
tively. Heavy lines correspond to m = 10~ 5 M Q yr _1 
and light lines to 1CT 6 M yr -1 . 

Fig. 4. — Properties of uniform clouds that are 
marginally stable to gravitational collapse in an envi- 
ronment with external pressure P = 2 x 1CP 10 dynes. 
Upper panel: Masses of neutral (light solid line) and 
molecular (heavy solid line) cores and including their 
corresponding HII envelopes (dashed lines) as func- 
tions of the critical core temperature T m . Lower 
panel: Density in core HI gas n (short dashed line), 
radii of HI core (solid line) and HII envelope (long 
dashed line) as functions of the critical core tempera- 
ture T m . 

Fig. 5. — X-ray spectrum at three radii within an 
HII- HI cloud having core mass Mhi = 3.5 Mq. The 
spectrum is shown at three radii: run = 0.206 pc 
(upper), tri = 0.0409 pc (middle), and r = pc 
(lower). 

Fig. 6. — Equilibrium structure of an HI core of mass 
Mhi — 3.5 Mq. Shown as functions of radius are 
the temperature T (long dashed line), HI density n 
(solid line), column density N (dotted line), degree of 
ionization x (dash-dot line), and the mass M (short 
dashed line). The pressure P — 8 x 10~ 10 dynes is 
uniform. 



Fig. 7. — Variation of the ratio of true mean temper- 
atures (T) within mass M(r) to the minimum temper- 
ature for gravitational support T m for five HI cores of 
different total mass. Regions for which (T)/T m < 1 
cannot support against gravity and collapse ensues. 
Each curve is labeled with the total HI core mass in 
solar units. 

Fig. 8. — Optical and X-ray surface brightness dis- 
tributions in NGC 4472: long dashed line: B-band 
surface brightness with photometric observations 
of Peletier et al. (1990) shown with open circles; 
short dashed line: ROSAT band surface brightness 
^rosat calculated from adopted density and temper- 
ature distributions; solid line: H/3 surface brightness 
E Hp (dropout) w ith q = 1 dropout model; dotted lines: 
H/3 surface brightness ^h/3(*) of gas ejected from stars 
for two assumed HII lifetimes, t* jC ond 

= 10 4 and 10 5 

years. 
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